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It has recently been found that Dirac and Weyl metals are characterized by an unusual weak-field 
longitudinal magnetoresistance: large, negative, and quadratic in the magnetic field. This has been 
shown to arise from chiral anomaly, i.e. nonconservation of the chiral charge in the presence of 
external electric and magnetic fields, oriented collinearly. In this paper we report on a theory of this 
effect in both Dirac and Weyl metals. We demonstrate that this phenomenon contains two important 
ingredients. One is the magnetic-field-induced coupling between the chiral and the total (or vector, 
in relativistic field theory terminology) charge densities. This arises from the Berry curvature and is 
present in principle whenever the Berry curvature is nonzero, i.e. is nonspecific to Dirac and Weyl 
metals. This coupling, however, leads to a large negative quadratic magnetoresistance only when the 
second ingredient is present, namely when the chiral charge density is a nearly conserved quantity 
with a long relaxation time. This property is specific to Dirac and Weyl metals and is realized only 
when the Fermi energy is close to Dirac or Weyl nodes, expressing an important low-energy property 
of these materials, emergent chiral symmetry. 


I. INTRODUCTION 

Weyl and closely related Dirac semimetals are the most 
recent addition to the growing family of materials with 
topologically-nontrivial electronic structure. Both were 
initially predicted theoretically, SH3 and realized very re¬ 
cently in a remarkable series of experiments.^®^ What 
distinguishes Weyl semimetals from other topologically 
nontrivial states of matter, like topological insulators 
(Ti) JHESDig that they are gapless. The topological object 
in this case is a point of contact between two nondegen¬ 
erate bands at the Fermi level, which acts as a monopole 
source of Berry curvature and thus carries an integer 
topological charge. The significance of such electronic 
structure featur es wa s emphasized in earlier pioneering 
work of VolovikpSB 7 which partly anticipated the recent 
developments. 

A hallmark of “topological” states of matter is the 
presence of metallic edge states, which arise necessarily 
due to the impossibility of a smooth connection between 
the topologically nontrivial sample and its topologically 
trivial environment. Weyl semimetal does indeed pos¬ 
sess such surface states, whose topologically nontrivial 
nature is manifest in the shape of their Fermi surface, 
having the form of an open arc (Fermi arc), rather than 
a closed curve, as in any regular two-dimensional (2D) 
metal. These have been seen directly using ARPES in 
the newly discovered Weyl semimetal materials TaAs and 

NbAsPira 

However, topologically nontrivial phases of matter of¬ 
ten also have unusual electromagnetic response, the most 
famous example being the precisely quantized transverse 
conductivity of a 2D quantum Hall liquid. Such a re¬ 
sponse is a robust, detail-independent manifestation of 
the nontrivial electronic structure topology on macro¬ 
scopic scales and is thus of particular interest. 

Topological electromagnetic response may be conve¬ 
niently expressed as a topological term, generated in the 


action of the electromagnetic field, when electrons in the 
occupied states are integrated out. In the case of the 2D 
quantum Hall liquid, this is the Chern-Simons term 

S = ~ [ dt<Pre vafi A v d a A f) , (1) 

47r J 

where (and henceforth) h = c = 1 units are used, and 
Filing factor one is taken above for concreteness. This 
may be generalized to 3D by stacking the 2D quantum 
Hall systems along a particular spatial direction,^ which 
gives 

S =~^2 j dtd 3 rG^^A v d a A 0 , (2) 

where G = 27r h/d is a reciprocal lattice vector in the 
stacking direction fi, corresponding to superlattice period 
d. Using integration by parts, Eq. (|2j) may be rewritten 
in the following form, which will prove useful below 

s = ^2 J dtd3r 0 ( r ) ( 3 ) 

where 9{r) = G • r. 

For a 3D TI, the topological term takes the form of the 
so-called 9- tenrPSI 

S = 3^2 fdtrfrOe^F^F af} , (4) 

where 9 = t: for a TI and 9 = 0 for a normal insulator 
(NI). Applying integration by parts as above, it is clear 
that Eq. Q in fact has no observable consequences in 
the bulk of the TI, since it is a total derivative and van¬ 
ishes upon integration by parts. It does have an effect 
on the TI boundaries: upon breaking time-reversal (TR) 
symmetry by e.g. magnetic impurity dopin g, it l eads to 
half-quantized anomalous Hall effect (AHE)f23E 25 ' 

Weyl semimetal may be regarded as an intermediate 
phase between a 3D TI (or NI) and the 3D quantum Hall 
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insulator, described by Eq. <[3|) .2] The topological term of 
a Weyl semimetal takes the form of Eq. ([3]), but with the 
field 0 given byPUlo] 

0(r, t) = 2b • r — 2& 0 t, (5) 

where 

b = b * = \Y. c ^- ( 6 ) 

i i 

Here C* is the topological charge of the ith Weyl node, 
K, is its location in momentum space and is its en¬ 
ergy. This is similar to topological terms arising in 
the context of Lorentz-invariance-violat ing extensions of 
the Sandard Model of particle physics .^ 1 As is easy to 
see, the linear space-time coordinate dependence of the 
0( r, t ) field is the only nontrivial dependence, compati¬ 
ble with space-time translational symmetry. Due to this 
space-time coordinate dependence, the topological term 
in Weyl semimetals does not vanish upon integration by 
parts, but instead takes the form, similar to Eq. <[2j) 

S = ~^2 j dtd 3 rd^ 9(r,t)e^A„d a A 0 . (7) 

This, in turn, leads to a nontrivial modification of the 
Maxwell equations in the bulk of the Weyl semimetal, 
which may be expressed as two extra contributions to the 
current density, obtained by varying Eq. ([TJ) with respect 
to the electromagnetic gauge potential 

3u = ^b^ va0 d a A p , M = 1,2, 3, ( 8 ) 

and 

> = -^boe 0va ^d a Ap. (9) 

Eq. ([ 8 ]) describes AHE with senri-quantized Hall conduc¬ 
tivity, proportional to the magnitude of the vector b, giv¬ 
ing the separation between the Weyl nodes in momentum 
space, while Eq. (0 describes the so-called chiral mag¬ 
netic effect (CMEjWl We will return to the meaning of 
the latter equation below. 

Topological term, describing the electromagnetic re¬ 
sponse of Weyl semimetals (and 3D TI as well), may be 
regarded as being a consequence of chiral anomaly ] 29 ! 33 ^ 37 
a fundamentally important concept in relativistic field 
theory, which has recently found its way into condensed 
matter physics and plays an important role in the mod¬ 
ern und erst anding of topologically-nontrivial phases of 
matter I2S1221 However, when using field theory concepts in 
the condensed matter context, one needs to exercise some 
care. Relativistic field theories possess exact symmetries, 
which in condensed matter systems may only be approx¬ 
imate. In particular, chiral anomaly is closely related to 
chiral symmetry, i.e. separate conservation of fermions 
of left and right chirality, which is an exact symmetry in 


theories of massless relativistic particles. Chiral anomaly 
refers to violation of this symmetry by quantum effects 
in the presence of electromagnetic field. In real Weyl or 
Dirac semimetals this symmetry may exist only approx¬ 
imately, if the Fermi energy is sufficiently close to the 
location of the nodes, so that the band dispersion may 
be taken to be linear to a good approximation. It is then 
unclear to what extent the concept of chiral anomaly is 
meaningful, when applied to Weyl and Dirac semimetals. 

In fact, the first sign of trouble with Eq. ([TJ) is its im¬ 
mediate consequence, Eq. 0, which describes the CME. 
This equation has the appearance of a current, driven by 
an applied magnetic field in the presence of an energy 
separation between the nodes. Such an energy separa¬ 
tion may exist in equilibrium in a noncentrosymmetric 
material,® which leads one to a problematic conclusion 
that Eq. (|9j) describes an equilibrium current. This would 
violate basic principles of condensed matter physics and 
can not happen . 41 The origin of this problem is precisely 
the relativistic invariance, assumed in the derivation of 
Eq. 0,® but not actually present in a real Weyl or Dirac 
semimetal. As was shown in Ref. l42[ in a condensed mat¬ 
ter setting, the response, described by Eq. depends 
on the order, in which the limits of zero frequency and 
zero wavevector are taken. When the frequency is taken 
to zero before the wavevector is taken to zero, which 
corresponds to thermodynamic equilibrium response, the 
current vanishes. However, when the order of limits is 
reversed, which corresponds to the DC limit of nonequi- 
librium response, the current is nonzero and given by 
Eq. (§. This dependence on the order of limits appears 
to violate Lorentz invariance and thus should not hap¬ 
pen in a relativistic particle physics context (at least if 
Lorentz invariance is assumed to be a fundamental sym¬ 
metry). This highlights the importance of being careful 
when using low-energy models, exhibiting “relativistic” 
properties, to describe Weyl semimetals and other Dirac 
materials. 

In this paper we will describe observable effects of the 
chiral anomaly, in particular the observable manifesta¬ 
tion of CME, in model Dirac and Weyl metals, i.e. lightly 
doped Dirac and Weyl semimetals. In accordance with 
the discussion above, we will use models for both, which 
explicitly do not possess chiral symmetry and are thus 
free of the artifacts of “relativistic” low energy models. 
We will demonstrate that in both Dirac and Weyl metals 
the main experimentally-observable consequence of CME 
is an unusual weak-field longitudinal magnetoresistance, 
which is negative, quadratic in the magnetic field, and 
large when the Fermi energy is sufficiently close to the 
Dirac or Weyl nodes.®l A shorter account of this work, 
devoted specifi cally t o the Weyl metal case, has already 
been published.®*^ The rest of the paper is organized 
as follows. In Section m we introduce the model we 
will use to describe both Dirac and Weyl metals, and 
which is based on the TI-NI heterostructure model of 
Weyl semimetals, introduced by us beforeP This model 
is the simplest model of a Dirac or Weyl metal that 
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does not suffer from the “relativistic” artifacts, in par¬ 
ticular it does not possess the spurious chiral symme¬ 
try of low-energy models with purely linear dispersion. 
In Section m we describe how CME manifests in Dirac 
metals. We derive coupled transport equations for, using 
field theory terminology, the vector and the axial (chi¬ 
ral) charge densities, which are coupled in the presence 
of an applied magnetic field. The coupling is shown to 
be the manifestation of CME. However, we demonstrate 
that this only leads to experimentally measurable conse¬ 
quences when a second ingredient is present: near conser¬ 
vation of the axial charge density, which is never exact, 
but becomes more and more precise as the Fermi energy 
approaches the Dirac node. In Section [TV] we describe 
the same effect in Weyl semimetals. The manifestation 
of CME in Weyl semimetals is found to be nearly iden¬ 
tical to the Dirac semimetals, i.e. whether the nodes are 
separated in momentum space or not does not matter 
for this effect. This appears to not be fully appreciated 
in the literature. We conclude in Section [V] with a brief 
discussion of our main results and experimental observ¬ 
ability of the effect. 


II. MODEL AND PRELIMINARIES 

We start from a model of Weyl and Dirac metals, based 
on TI-NI multilayer heterostructure, introduced by usP 
The advantage of this model is that it is extremely sim¬ 
ple, yet more realistic than the most generic low-energy 
model of a Dirac or Weyl metals would be, in particular 
it does not have the unphysical chiral symmetry. Since 
the model has already been described in a number of 
publications, here we will only recap the most essential 
points. The momentum space Hamiltonian, describing 
the multilayer structure, is given by 


breaking term ba z , which arises physically either from 
polarized magnetic impurities or an external mag netic 
field. This has been described in our earlier papers,®^ 
and we will not dwell on it further. 

To make contact with chiral anomaly, it is useful to 
recast Eq. (10) in a “relativistic” form. To this end we 
expand Eq. (10) to leading order in the crystal momen¬ 
tum near the point k = ( 0 , 0 , 7 r/d), at which the gap 
closing occurs at the TI-NI transition when A s = An- 
We obtain 


H = vft z (z x er) • k + A ndr v k z + (As — An) t x . (12) 

This implies the following representation of the first four 
Dirac gamma matrices 

0 x 1 • y y 2 • y x 3 ■ z f-i 

7 = r , 7 = zr y a y , 7 = —iT y a ,7 = it . (13) 


The fifth gamma matrix 


q 5 = = t v (J z 


(14) 


defines the axial charge operator n a = 7 5 = T v cr z . Ab¬ 
sorbing the Fermi velocities Vp and vp = And into the 
definition of the corresponding momentum components 
and replacing k y -A —id^, one obtains the following “rel¬ 
ativistic” Lagrangian 


C = ip^id t ip - H = ip + ieA^ + ib ^q 5 ) - to] ip. 

(15) 

Here to = As — Ad is the Dirac “mass”, ip = ^q° is 
the Dirac adjoint of the Grassmann field ip, A „ is the 
electromagnetic gauge potential, and we have also intro¬ 
duced chiral gauge field b M . Explicitly, chiral gauge field 
arises from the following terms, added to the Hamilto¬ 
nian Eq. (12 1 


H b = b 0 r y a z + b\T x (r x + b 2 r x a y + b 3 a z . 


(16) 


H = vft z (z x er) • k + A (k z ). (10) 

Here z is the growth direction of the heterostructure, 
vf is the Fermi velocity, associated with the motion in 
the transverse (x, y) directions, er are Pauli matrices, de¬ 
scribing the real spin degree of freedom, while r is the 
pseudospin, describing the top and bottom surfaces of 
TI layers in the heterostructure. The operator A (k z ) de¬ 
scribes the electron dynamics in the growth direction and 
is explicitly given by 


A (k z ) = A s t x + ^ (r + e ikzd + h.c .) , (11) 

where A s,n are amplitudes for tunnelling between top 
and bottom surfaces of the same (S) or neighbouring (D) 
TI layers and d is the superlattice period. We will take 
them both to be positive for concreteness. This structure 
is an ordinary insulator when A$ > A n and a strong 3D 
TI otherwise. The point As = A n marks the TI-NI 
phase transition. At this point the structure is a Dirac 
semimetal. Weyl semimetal is obtained by adding a TR 


The first term in Eq. (16) is clearly an axial chemical po¬ 
tential term, which shifts the left (L) and right-handed 
(R) components of the Dirac fermion in opposite direc¬ 
tions in energy. The last term is magnetization (or mag¬ 
netic field) in the ^-direction, which shifts the L and R 
components in opposite directions in momentum space 
along the 2 -axis. The second and third terms have the 
same symmetry as magnetization components in the x , y 
directions, and thus may be regarded as such. However, 
one needs to be aware that bare a x,v operators will have 
a very different effect on the spectrum, creating a nodal 
line state rather than point nodes 

Chiral anomaly refers to anomalous nonconservation 
of the axial current Jpp = ip^ IJ -'y 5 ip. This means that the 
axial current continuity equation has the following form 


du. = 2 im + ^e^F^F af3 . (17) 

The first term on the right hand side of Eq. ( fl7| ) is the 
classical contribution to the axial charge continuity equa¬ 
tion, which is easily obtained from the Dirac equation. 
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When the mass m = 0, i.e. when As = Ajj, this term 
vanishes. This is an expression of chiral symmetry, i.e. 
classically the axial charge is conserved when m = 0 . 
This conservation is violated when the continuity equa¬ 
tion is evaluated in the second-quantized theory, which 
is where the second term comes from. However, m = 0 
when As = Ap is only obtained at leading order in the 
expansion of the operator A (k z ) near k z = n/d. In fact, 
to is a function of k z and does not generally vanish. This 
means that the chiral symmetry is always only approx¬ 
imate and the axial charge in never a truly conserved 
quantity, even when the anomaly term is neglected. It is 
then clear that anomaly-related effects may only be ob¬ 
servable if chiral symmetry is almost there, i.e. the axial 
charge relaxation time is long. We generally expect it to 
be long when the Fermi energy is close to zero, and thus 
the effect of higher order terms in the expansion of A(fc~) 
is not significant. The purpose of the rest of the paper is 
to make these statements quantitative and evaluate the 
axial relaxation time explicitly for model Dirac and Weyl 
metals, described by Eq. (10 1 . 


III. ANOMALOUS DENSITY RESPONSE IN A 
DIRAC METAL 


In this section, starting from a microscopic model of a 
Dirac semimetal, given by Eq. (10), we will derive trans¬ 


port equations for the vector and axial charge densities, 
in the presence of an external magnetic field. As will 
be demonstrated, chiral anomaly will manifest in these 
equations as a coupling between the vector and the ax¬ 
ial charge densities, induced by the magnetic field. We 
will discuss under what conditions this coupling leads to 
observable transport phenomena, namely quadratic neg¬ 
ative magnetoresistance. 


We start from the Hamiltonian Eq. (10), with an added 
magnetic field in the z-direction 


H = vpT z (z x er) • (— i'V + eA) + A (k z ). 


(18) 


index a is a composite index a = s,t, where t = ± la¬ 
bels two components of the Kramers doublet. The en¬ 
ergy eigenvalues do not depend on t as we have ignored 
the Zeeman splitting due to the applied field. tA{k z ) 
are the two eigenvalues of the A (k z ) operator, where 
A{k z ) = y/A% + A 2 d + 2A s A d cos (k z d). The corre¬ 
sponding eigenstates have the following form 

\n,a,k y ,k z ) = Y [z^ T (k z )\n- l,k y ,k z ,t,r) 

T 

+ Z nlr( k z)\n,ky,k z ,l,T)] , (21) 

where 

(: r\n,k y ,k z ,a,T) = -^=e lk * z (t> n k y (r)|cr, r), ( 22 ) 

4> n k y (r) are the Landau-gauge orbital wavefunctions and 
< 7 , r are the spin and the top-bottom surface pseu¬ 
dospin labels respectively. The four-component eigen¬ 
vector | z“(k z )) may be written as a tensor product of 
the two-component spin and pseudospin eigenvectors, i.e. 
\ z n( k z)) = |<(&*)) ® | u a {k z )), where 


K\k z )) = 

!«*(**)> = 


1 

71 

i 

71 



tA(k z ) 

€n(k z )’ 


—is 



A s + A D e~ ikzd 
A (k z ) 


1 jA(k z ) \ 
e n (k z ) ) ’ 


(23) 


The lowest n = 0 LL is special, which is a consequence 
of nontrivial Berry curvature. The s quantum number is 
absent in this case and taking B > 0 for concreteness, we 
have 


£nt{k z ) = -tA(k z ), (24) 

and | vl(k z )) = ( 0 , 1 ). 

We model the potential due to random impurities as 
Gaussian white noise potential with (V(r)) = 0 and 


We will adopt Landau gauge for the vector potential 
A = xBy and ignore the Zeeman splitting in this section. 
The restriction of the magnetic field to the ^-direction 
simplifies calculations in the context of our model, but is 
otherwise nonessential. After a canonical transformation 


a z 


(19) 


Eq. (18 1 is easily diagonalized. The eigenvalues have the 
form 

ena(k z ) = s^2uj 2 B n + A 2 {k z ) = se n (k z ). (20) 


Here n > 1 is the main Landau level (LL) index, k y is 
the intra-LL orbital label, s = ± labels the two sets of 
positive and negative energy (or electron-like and hole¬ 
like) eigenvalues, ujb = Vf/^b is the Dirac cyclotron fre¬ 
quency and £b = 1 /VeB is the magnetic length. The 


(V’(r)V'(r')) = 7 2 <5(r — r'). (25) 

For simplicity we assume that the impurity potential is 
independent of the pseudospin index r, which physically 
means that we ignore scattering between the top and 
botton surfaces of the TI layers. This does not affect the 
results qualitatively and is only done to simplify calcula¬ 
tions. 

The impurity scattering may be treated by the stan¬ 
dard diagrammatic perturbation theory. In the self- 
consistent Born approximation (SCBA) one obtains the 
following expression for the retarded impurity self energy 
(see Fig. [l] (a) for graphical representation) 

(k z ,u) = ^~ {\{nakyk z \V\n'a'k' y k' z )\ 2 ) 

2 n'a'k' y k' z 

x G* 0 ,(fc», 


(26) 
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where 


G*(*„w) = 


W - £na(k z ) + if]' 


(27) 


is the retarded Green’s function of a clean Dirac metal 
(full self-consistency is unnecessary here) and £, na (k z ) = 
e na (k z ) — e F . We will assume that the Fermi energy 
e F > 0, i.e. that the Dirac semimetal is electron-doped. 
Matrix elements of the impurity potential between the 
LL eigenstates are easily evaluated and are given by 

(nakyk z \V\n'a'k' y k' z ) = ] V 

X V(q)^, fea _ fe ^ 9zifcz _ fei e^^^+ fc v)/ 2 
x 'y 1 ^ z n J fT(kz)z n i^ T (k z )F n -i in i-i(q) 

T 

+ Zn iT {k z )z^ iT {k' z )F n ^( q) , (28) 

where the LL formfactors have the following well-known 
form 


n . fn^ f iqx^B ~ q y i B 

F "’"' (q) v nT V- 2 - 


X e- q2e °/ A L n n 7 n ' ( 


n 2 f 2 

q *-B 


(29) 


and Z£r" (x) are the generalized Laguerre polynomials. 
Using the following properties of the LL formfactor mo¬ 
mentum integrals 


( 2 7 r ) 2 ii, "."' (q)F "'."( q) - 2 7 rf|’ 


d 2 q 

Wf‘ 


Fn,n' ( q)Tn'_ l,n— 1 ( q) — 0; 


(30) 


we obtain 

r 

= 7 


na,n' a‘ 
2 


(k z ,k' z ) = (\(nak y k z \V\n'a'k' y k' z ) | 2 ) 

2 | ..a 7 


<t(^)IK't(^)l' 2 + K|(fc,)rKq(^)h 

X \(u a (k z )\u a '(k ' z ))\ 2 . 

The SCBA equation then takes the form 


(31) 


zR 


(k z ,U)) — L y ' ^na,n'a'(kz,k z )G n , al (k z ,u). 

B z * ■ ’ * 


( 32 ) 

At this point we will assume that e F is sufficiently large, 
so that scattering between the electron- and hole-like LL 
may be neglected. Then the negative energy s = — states 
do not contribute and we will ignore them henceforth. 
We will also drop the explicit s = + index, since all 
states are the s = + states, and replace a index by t 
from now on. 

The SCBA equation may now be solved analytically. 
Using 


JmG„ t (k z ,u) = -nS[e n (k z ) - e F ], 


(33) 


+— = —>— + —*—£—:—► 

(a) 


-*-1 


r+- ► 


r¥- 



= ■ + ■ 





U- -4-M- < 4 


b*- 


(b) 


FIG. 1. (a) Graphical representation of the SCBA equation. 
Regular lines denote bare Green’s functions, while bold lines 
correspond to impurity-averaged Green’s functions. Dashed 
line denotes impurity averaging, (b) Graphical representation 
of the diffusion propagator, which is the sum of all ladder di¬ 
agrams. Bold lines correspond to impurity-averaged Green’s 
functions. 


it is easy to see that the dependence of the matrix ele¬ 
ment T n t, n 't> {k z , k' z ) on the LL indices n , n' in fact drops 
out, since this dependence only enters through the LL 
energies e n (k z ), which may simply be replaced by the 
Fermi energy. We then obtain the following expression 
for the impurity scattering rate 


- 77 TU = - 2 ImS nt( fc ^ w ) 
T \ K z) 

A s + A D (cos {k z d)Y 
e F 


1 

TO 


1 As + Ap cos (k z d) 
e F 


(34) 


Here 


— = xtt7 2 g{e F ), 

To 2 


(35) 


and g(e F ) is the total density of states at Fermi energy. 
Explicitly 


9^f) = 


2t tC 


B J—n/d 


n / d rib . 

-^^2 s [ e n(kz) ~ e F )] 


e F 


T / d dk 


= -V/ wf©[^-A (k z )}, 


' —n/d 


27r 


(36) 


where in the second line above we have assumed that the 
magnetic field is weak and converted the sum over the 
LL index n to an integral. Finally, (cos {k z d)) in Eq. (34) 
means the average of cos (k z d) over the Fermi surface, 
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which is defined as 


o r n / d r iu w _, 

(cos(M)> = —— / V cos (k z d)S[e n (k z )-e F )\. 

9 \ e F) J—n/d 27 r ^ 

(37) 

We will use this definition of Fermi surface averages 
throughout the paper. Evaluating the average in Eq. (37) 
in the weak-field limit, one obtains 


(cos (k z d)) 


1 /l ( + ^ d e F \ 
MV \ 2 AsAd ) 


(38) 


where 


ko = - arccos 
d 


( £ f \ 

\ 2A s A d ) ’ 


(39) 


is the solution of the equation A (k z ) = 6f modulo ir/d. 

At this point we will specialize to the case of a Dirac 
metal, i.e. set A s = Ad, which makes the Dirac mass 
at the point k x = k y = 0, k z = tt/ d vanish. We will 
also assume that the Fermi energy is close to the Dirac 
point, which means ep/ Ag <C 1 (but still far enough that 
cft o 1). Then we obtain 


(cos (M)) 



180A| 




(40) 


and it is clear that, to leading order in the small param¬ 
eter €p/As, 


1 _ 1 
r(k z ) t 0 ' 


(41) 


We now want to find the propagator of the diffusion 
modes in our system, which correspond to nearly con¬ 
served physical quantities with long (i.e. much longer 
than To) relaxation times. In the limit epTo 1, the dif¬ 
fusion propagator may be evaluated by summing ladder 
impurity scattering diagrams, as shown in Fig. [l] (b)P^I 
This approximation (self-consistent non-crossing approx¬ 
imation) is consistent with the SCBA for the impurity 
self-energy, in the sense that together they preserve ex¬ 
act conservation laws, in particular charge conservation. 
The diffusion propagator (or diffuson), evaluated in the 
self-consistent non-crossing approximation, takes the fol¬ 
lowing form 


D(q,D) = [l-/(q,D)]- 1 , (42) 

where is a 16 x 16 matrix with respect to the 

combined spin and pseudospin indices a = cr,T, which 
has the following explicit form 

2 p 

Ia iaa ,a 3 a 4 (q,fi) = / d\d\' 

xG« a3 (r,r'|D)G^ 4Ct2 (r',r|0). (43) 


The Green’s functions, appearing in Eq. (43) are the 
impurity-averaged SCBA retarded and advanced Green’s 
functions, found above 


G^( r ,r'|Q)= Y, 

ntk,,k z 


where 


(ra | ntky k z ) ( ntk y k z 


| r'a') 


kl - £, nt { k z ) ± ?’/2to 


(44) 


(ra\ntk y k z ) = -j=e lk * z z^k^n^r). (45) 

Evaluating /(q, D) in the general case is a formidable 
task, mostly due to the presence of the LL index sums. 
LLs with different n will generally be mixed by impu¬ 
rity scattering, in which case analytical evaluation of 
Eq. (43) becomes impossible. However, we are primarily 


interested in transport along the direction of the applied 
magnetic field, i.e. the ^-direction. If, in accordance with 
this, we set q = qz , it is easy to see that only a single LL 
index sum remains in Eq. (43) after integration over the 


x, y coordinates. In this case we obtain 


L O'! OL2 , 0 : 30:4 


(<?,fi) = 


7 

2-kP b L 


E 

ntt'k z 


z t nai {kz)z t na3 {k z ) 

G — £, nt .{ k z + q/2 ) + */2 tq 


A 04 (k z )z n ( k z ) 


~£,nt’{k z - q/2) - i/2r 0 ‘ 

The dependence of the eigenfunctions z^ lct (k z ) on q has 
been neglected in Eq. (461, since the corresponding terms 


(46) 


are subdominant in the limit epT 0 1 and qvpT 0 <C 1. 

At this point we need to explicitly separate out the part 
of the diffusion propagator that corresponds to hydro- 
dynamic modes, i.e. modes with long relaxation times. 
On physical grounds, we expect only two such modes 
to be present in our system, corresponding to the diffu¬ 
sion of the vector n v = a°T 0 and the axial n a = a z r v 
charges. Note that the axial charge operator changes to 
n a = <j°t v after the canonical transformation Eq. (19). 


The projection onto the vector-axial charge subspace is 
accomplished by the following transformation 


Kb 1 


= -(err' 
4 v ' 


V 


-1 


CX. 20 t\ OL\CX -2 , 0 : 30:4 


(err 


) b 

' 0 : 30:4 ’ 


(47) 


where a, b refer to either the vector or the axial charge 
and summation over repeated indices is implied. The 
projected inverse diffusion propagator is a 2 x 2 matrix. 
Its diagonal components describe the independent trans¬ 
port of the vector and axial charge densities, while the 
off-diagonal components describe their coupling, induced 
by the applied magnetic field. Let us first evaluate the 
off-diagonal component (the two off-diagonal components 
are equal by reciprocity). We have 

= -\( <jQT °)cx2 ai I ai a 2 ,a 3ai { c l- ^)( 0 -V ) a3Q4 

= 7 2 (4'(MI4(fcz)) 

87 t£ 2 b L z Q - ^ nt (k z + q/2) + i/2r 0 

TlsL'L K z 

x (Zn(k z )\T y \z^(k z )) 

~£,nt' (k z — q/2) — i/2r 0 ' 


(48) 
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Since {z t n (k z )\z t n (k z )) = 5 tt >, and 
{z\k z )\Ty\z t {k z )) = -t 


Ad sii \{k z d) 
A (k z ) ’ 


(49) 


it is clear that LLs with n > 1 can not contribute to 


Eq. (481, since their energies e nt (k z ) do not depend on 


the index f, which leads to an exact cancellation of con¬ 
tributions with t = ±. The n = 0 LL, on the other hand, 
does contribute, since the corresponding eigenstate ener¬ 
gies do depend on t, as seen in Eq. (24). When e F > 0, 


it is clear that only the t = — lowest LL contributes to 
Eq. (48), since only this LL crosses the Fermi energy. 


Using an identity 


AB = 


B — A 


A -1 — B~ x ’ 


(50) 


we have 


1 


1 


x 


n - Co -(kz + q/2 ) + i/2r 0 -Co ~{k z - q/2) - i/2r 0 
1 

n - Co -(kz + q/2 .) + Co— {kz - q/2) + i/r 0 
1 1 
-Co -(kz - q/2) - i/2r 0 fl - Co -ik z + q/2) +i/2 r 0 _ 


2ttz 5[A(k z ) — €f\ 
n + *Ao 


Substituting this into Eq. (48) and expanding to first 
order in 12 and q , we obtain 


z\,«(q,= 


7 T o 


pn / d 


dky 


1 dA 


87t£| z A s d dk 


d[A(fc,) - e F ] 


I r, dA 
x ( 1 + zIItq - ig—r 0 


(52) 


where we have used 
dA AsAp,dsm(k z d) 

dk z A{k z ) 


=-A s d(z {k z )\r v \z ( k z )). 

(53) 

Since dA /dk z is an odd function of k z with respect to the 
point k z = 7r/d, it is clear that only the term, propor¬ 
tional to q , survives the integration over k z in Eq. (52). 
Then we obtain 


v V a{ q, = -k 


7 2t o 


4.n£%A sd 


dA 


dk. 


(54) 


k z =kt 


where kj = it/d±ko are the two solutions of the equation 
A(k z ) = e F . Evaluating the lowest LL Fermi velocity 
dA/dk z explicitly, we obtain 


dA 


dk. 


k,=kt 


= -U4A 2 - 


A s d. 


(55) 


Then, using y 2 = 2/Trg(e F )T Ql we finally obtain 


V V a{ q, 12) = 12av (q, °) = -iq T o 


eB 

2TT 2 g(e F ) 


= -iqT 0 T. 


(56) 


The coefficient T in Eq. (561 is a new transport coef¬ 


ficient, that describes the lowest LL-mediated coupling 
between the vector and the axial charge densities. This 
coupling may be regarded as being a consequence of chi¬ 
ral anomaly. 

The diagonal elements of the inverse diffusion prop¬ 
agator correspond to independent transport and relax¬ 
ation of the vector and axial charge densities. These are 
nonzero in the absence of the magnetic field and we will 
thus evaluate them in the limit B —» 0, since we are in¬ 
terested in weak-field transport here. Accordingly, the 
contribution of the lowest LL is negligible in this case 
and we will ignore it. The limit of B —> 0 will be taken 
after summing the contributions of all the n > 1 LLs. 
One obtains 


**(*«>-i- 01 : r 

„>i J-*/d 


dk z 


S[e n {k z ) - e F \ 
1-iClro + iqTo^j^' 

, (57) 

The sum over the LL index n may be done in the limit 
B —> 0 by converting the sum into an integral, just as was 
done when solving the SCBA equation above. Perform¬ 
ing the integral and expanding to leading non vanishing 
order in 12 and g, one obtains 


(51) 




(58) 


where the angular brackets denote average over the Fermi 
surface, defined as in Eq. (37). Using 


A dA IdA 2 . . . . 

A *: = 2*7 = “ AsADSm<M) ’ 

the average is easily evaluated and we obtain 

((«)»■«• 


(59) 


(60) 


Defining the 2 -direction diffusion coefficient as D = 
( Asd) 2 r 0 /3 , we finally obtain 


1> v v( q, ty = -212t 0 + Dq-T 0 . 


(61) 


This has the expected form for the inverse diffusion prop¬ 
agator of a conserved quantity. Namely, the full diffu¬ 
sion propagator 2?(q, 12 ) will exhibit a diffusion pole at 
12 , q —)■ 0 as a consequence of an exact conservation of the 
vector charge. 

Finally, we need to evaluate V~2. Here we expect that 
the diffusion pole will be absent due to a finite relaxation 
rate for the axial charge density, since it is not an exactly 
conserved quantity. We obtain 


V~2(qM) = l- 


7 2 t° 


E 


f*7T/d 


87r£p — J-Tr/d, 
B n>l,tt' 7r /“ 

S[e n (k z ) - e F \ 

l-mTo + iqTo-^^’ 


dfc 2 |(4(MIU/l4 (k z ))\ 2 


( 62 ) 






































Evaluating the matrix element in Eq. 


one obtains 
2 / 


(kJ) . 

1 w e ^ Kz > 

(63) 

Substituting this back into Eq. (62) and evaluating the 


sum over the LL index by converting it to an integral, as 
before, and expanding to leading non vanishing order in 
and q, we get 

Vaa( q,fi) = 


1 - ( 1 - 


A^-A^sin 2 (M)\ (1+ . nro) 


■ 2 2 
+ Q r o 


1 - 


A 2 (k z ) — A 2 D sin 2 (fc 2 c?) ^ ^ A dA 

€f dk z 


(64) 


where the angular brackets again mean average over the 
Fermi surface. Evaluating the Fermi surface averages, 
assuming as before that cf/As <C 1, we finally obtain 


£>aa(q> = -i^To + — + Dq 2 r 0 , (65) 

Ta 

where 



is the axial charge relaxation rate. Eq. ( |66| is one of 
the main results of this section. As expected, the axial 
charge is not exactly conserved, as the chiral symmetry 
is always explicitly violated in a real Dirac semimetal by 
nonlinearity of the band dispersion, which is necessar¬ 
ily present. However, since the band dispersion becomes 
more and more linear as the energy is reduced towards 
the Dirac point, the axial relaxation rate tends to zero 
as the Fermi energy goes to zero faster than the momen¬ 
tum relaxation rate 1 / to , which of course also vanishes 
in the limit cf —> 0 due to the vanishing density of states 
(This is true provided we neglect the influence of the 
magnetic held on the density of states. In principle, even 
in the limit €f —> 0 there is a finite density of states, pro¬ 
portional to B. We ignore this in the weak-Held limit). 
Thus, near the Dirac point r a t 0 , which expresses the 
near-conservation of the axial charge due to the emergent 
low-energy chiral symmetry. As will be seen below, this 
is a necessary condition for a large negative magnetore¬ 
sistance. 

Collecting all the matrix elements, we obtain the fol¬ 
lowing result for the full inverse diffusion propagator, 
which describes coupled transport of the vector and axial 
charge densities 


performing the inverse Fourier transform, we obtain the 
coupled diffusion equations 


dn v _ d 2 n v dn a 

dt dz 2 dz 

dn a _ d 2 ?z a _ ria „ dn v 

dt dz 2 r a dz 


( 68 ) 


Since the vector charge is exactly conserved, the first 
equation must have the form of the continuity equation 


dn v 

~dt 


^ ' jlM 


(69) 


where j„ is the vector current, i.e. current of the vector 
charge. This leads to the following explicit expression for 
the vector current 


jv — 


Q-q d/i v 
e dz 



(70) 


Here and /x a are the vector and axial electrochemical 
potentials, er 0 = e 2 g(e^).D is the zero-field Drude con¬ 
ductivity, and we have used <5n u>a = g(tF)^ v ,a- The 


first term in Eq. (70) is an ordinary current in response 
to a gradient of the electrochemical potential. The sec¬ 
ond term is a consequence of chiral anomaly and is an 
extra contribution to the current, proportional to the 
applied magnetic held and (nonequilibrium part of) the 
axial electrochemical potential. This is known as CME 
in the literature,® and this extra contribution is what 
leads to the anomalous negative longitudinal magnetore¬ 
sistance. However, CME by itself is only one component 
of the experimentally observable effect, i.e. the negative 
magnetoresistance. The second crucial component, with¬ 
out which the effect is unobservable, is contained in the 


second of Eq. (68). Namely, as will be seen shortly, the 


CME leads to observable magnetoresistance only if the 
axial charge relaxation rate l/r 0 is small, i.e. the ax¬ 
ial charge is a nearly conserved quantity. As discussed 
above, this near-conservation of the axial charge is a char¬ 
acteristic feature of Dirac (and Weyl) semimetals, which 
becomes more and more precise as the Fermi energy is 
reduced towards the Dirac (or Weyl) point. 

To obtain the CME-related magnetoresistance explic¬ 
itly, we now assume that there is a uniform steady state 
vector current density in the z-direction j v , present in 
the system. The second of Eqs. (68) then gives 


, -p d/i v 

dfj, a = r T a ——. 

oz 


(71) 


Substituting this into the equation for the vector current 
Eq. (70), we obtain the following expression for the total 


diagonal conductivity 


v~ i (q,n) = 


—i£lTQ + Dq 2 TQ —iqTro 

-iqTT 0 -if2r 0 + r 0 /r a + Dq 2 r 0 


(67) 


Viewing Eq. (67) as the inverse Green’s function of the 
diffusion equation for the vector and axial charges and 


c zz 


e 4 B 2 T a 

a ° + 4t T 4 g(e F )' 


(72) 


Thus CME is manifested as a positive longitudinal mag¬ 
netoconductivity (or negative magnetoresistivity), pro¬ 
portional to B 2 . Crucially, it is also proportional to 
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T a , and a large r a is thus necessary for this effect to be 
significant. The magnetoresistance is further enhanced 
when e F —> 0 by vanishing of the density of states as 
g(e F ) ~ e 2 F . 


IV. ANOMALOUS DENSITY RESPONSE IN A 
WEYL METAL 


In this section we will extend the theory of the 
anomaly-related negative magnetoresistance, presented 
above, to the case of Weyl metals, where the individ¬ 
ual Weyl fermion components of the Dirac fermion are 
separated to distinct points in momentum space. A 
shorter account of this work has already been presented 
in Refs. 04] and 05] As the calculations are quite simi¬ 
lar to the case of Dirac metals, described in the previous 
section, here we will only focus on the differences from 
the Dirac metal case and skip some of the details. 

In the context of our model Dirac semimetal, described 
by Eq. (18), the separation of the Dirac fermion into 


Weyl fermions is most easily accomplished by adding a 
term bcr z to the Hamiltonian. Physically this term may 
arise from magnetized impurities, doped into the Dirac 
semimetal, or even from the Zeeman coupling to the ap¬ 
plied magnetic field. 

The LL energy eigenvalues now have the form 


i{k z ) = s\J2u B n + m 2 {k z ) = se nt {k z ), n > 1, (73) 


while 


eo t(k z ) = —m t (k z ). 


(74) 


Here 


m t {k z ) = b + tA(k z ). (75) 

Taking b to be nonnegative, m-(k z ) vanishes at two 
points along the 2-axis in momentum space, given by 
the two solutions of the equation 


A (k z ) = b. 


The solutions are k z = ir/d ± fco, where 


ko = - arccos 
a 


/ A| + A 2 D -b 2 \ 
V 2,A F A D ) 


(76) 


(77) 


These correspond to the locations of the two Weyl nodes 
of opposite chirality on the 2-axis in momentum space. 
The nodes exist as long as b c \ < b < b c2 where b c i = 
| As — As>| and b c2 = As + Ap,. The eigenvectors are 
given by 


\<(k z )) = 
I u\k z )) = 



m. t {k z ) 


£nt{kz) 
A D e~ ikzd 

A(kJ 


, —is.\ 1 — s 


m t (k z ) 
€nt{kz) 


(78) 


while the n = 0 LL is polarized downwards, as before. 

The main difference from the Dirac metal case is that 
the Kramers degeneracy between the t = ± states is now 
broken by the spin splitting term ba z . When b is suf¬ 
ficiently large (i.e. b > e F ), we may ignore the t = + 
states entirely. Solving the SCBA equations as before, 
we obtain 


1 

T(k z ) 


1 

TO 


1 + 


to (k z )(m-) 


where 1 /tq = nj 2 g(e F ) and 


9^f) 


1 

2n£ 2 B 



dk z 

27T 


^2S[e n -(k z ) - e F ], 


(79) 


(80) 


is the total density of states at Fermi energy. At the 
Weyl nodes m-{k z ) vanishes and changes sign. This im¬ 
plies that, for sufficiently small Fermi energy, such that 
the band dispersion in the 2-direction may be assumed 
to be linear to a good approximation, the Fermi surface 
average (to_(/c.)) will vanish. This property has a sim¬ 
ple geometrical interpretation. The Weyl nodes may be 
thought of as monopole sources of Berry curvature, whose 
2-component is proportional to m-(k z ). This clearly av¬ 
erages to zero when integrated over the volume, enclosed 
by a sufficiently small Fermi surface sheet, containing the 
nodeP^I Assuming this to be the case, we obtain 


1 ~ 1 

r{k z ) t 0 ' 


(81) 


Note that since the densities of states in Eqs. (36) and 
(80) are essentially identical in the limit of small Fermi 
energy (the two-fold Kramers degeneracy in the Dirac 
metal case is replaced by two identical Fermi surfaces, en¬ 
closing the Weyl nodes, in the Weyl metal case), the im¬ 
purity scattering rate in the Weyl metal is twice as large. 
This is easy to understand physically and follows simply 
from the near orthogonality of the |u ± (fc 2 )) eigenstates 


at small momentum difference, i.e. (n*(fc 2 )|M* (fc 2 )) ~ Stt' 
when |fc 2 — k' z \d 1, which means that scattering be¬ 
tween the two components of the Kramers doublet is sup¬ 
pressed in the Dirac metal. This suppression disappears 
in the Weyl metal case, since in this case the two Fermi 
surfaces arise from states in the same t = — band. 

The evaluation of the diffusion propagator goes along 
exactly the same lines as before. The only difference is 
that only the t = — states contribute, i.e. we have 


'010:2,0:304 


(«,«) = 




"noi (.kz^ri 


1 (kz) 


2n£ 2 B L z ^ fl - in-{k z + q/2) + i/2r 0 


Z nai ( k z) Z n a2 (kz) 
-£, n -(k z - q/2) - i/2r 0 ' 


(82) 


For the same reason, the projection of the full diffusion 
propagator onto the vector and axial charge subspace 
differs by a factor of 1/2 from the Dirac semimetal case 


V 


-1 

ab 


1 

2 


(err) 


a 

&2&1 


V 


-1 

Q' 1 a 2 , a 3 a 4 


(err) 


b 

0:30:4 7 


(83) 
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The vector to axial charge coupling term arises, as before, 
entirely from the contribution of the n = 0 LL. We obtain 


Vvai q> 


2nf%g(e F )A s d 


dA 

dk z 


k z =k± 


(84) 


where kf are now solutions of the equation 


A(fc z ) = b + e F , (85) 


which determines the points at which the Fermi energy 
intersects the n = 0, t = — LL. One obtains 


= 0(h t J [( b + ^) 2 - b ciWc2-(b + e F n 

k z =ki 2(6 + e F )V 

( 86 ) 

Assuming b c \ <C b + e F <C b c2 , which implies that the 
Weyl node splitting and the Fermi energy are such that 
the band dispersion at the Fermi level may be taken to 
be linear, one obtains 


dA 

dk z 


dA 

dk z 


k z =k : t 


As + Ap 
2 


d 


A s d. 


(87) 


This gives 


V va (q. 2^g( y e F ) ~ ~ iqT ° T ’ ( 88 ) 

i.e. an identical result to what we obtained before in the 
case of the Dirac metal. 

The diagonal elements of the inverse diffusion propa¬ 
gator are also evaluated in exactly the same way as in 
the case of the Dirac metal. The form of the expression 
for the vector charge part of the propagator is, as before, 
constrained by the vector charge conservation 


i.e. again identical to the corresponding Dirac metal re¬ 
sult, obtained in Section [TTT1 

Finally, for the axial charge block of the inverse diffu¬ 
sion propagator we obtain the following expression 

Vaa (q- V) = -ittT 0 + — + Dq 2 T 0 , (93) 

7"a 


where the axial charge relaxation rate is now given by 


to _ 1 - (v f /A s d) 2 
r a ( v F /A s d) 2 


(94) 


This expression for the axial charge relaxation rate rep¬ 
resents the most significant difference of the Weyl metal 
case from the Dirac metal case. This equation shows, 
in particular, that in a Weyl metal the (dimensionless) 
axial charge relaxation rate is essentially always finite, 
even in the limit e F — > 0. It may still be expected to 
be small, which is easily seen explicitly in the limit when 
b c i b <C b c 2, and e F -C b. In this case we obtain 


1 

T a 



(95) 


i.e. the axial charge relaxation increases quadratically 
with the spin-splitting parameter b as the Weyl nodes get 
split and separated out of the parent Dirac metal state. 
The diffusion equations themselves and the magnetoresis- 
tance formula are identical to the Dirac semimetal case, 
with the axial charge relaxation rate given by Eq. (951, 
so we will not repeat them explicitly here. 

Before we conclude this section, we would like to point 
out an important caveat, which applies to the results of 
this section. Namely, all of the results above are only 
applicable if the condition 


b c i < b < b c2 , 


(96) 


Zw(q, O) = -iDr 0 + Dq 2 r 0 , 


where the diffusion coefficient is given by 


D = v f t 0 


i-(k z ) 


(89) 


(90) 


which is identical to the corresponding result in the Dirac 
metal case but with A (k z ) replaced by m_(/c z ). The 
Fermi velocity in Eq. (90) is v F (k z ) = \dA/dk z \ } evalu¬ 
ated at the Weyl node locations, which is given by 


v F = 


2b 


\J(b 2 -b 2 cl )(b 2 c2 -b 2 ). (91) 


The average of m_(k z ) over the Fermi surface may be 
easily evaluated in the limit of small Fermi energy, which 
allows one to expand m_(fc z ) to leading order in devia¬ 
tion of k z from the Weyl node locations. In this case one 
obtains 


D « ^v 2 f t 0 , 


(92) 


is satisfied. An implicit assumption here is that Weyl 
semimetal is obtained from a parent state, which is nearly 
a Dirac semimetal (hence the small b c i = |As—Ao|), and 
we are far from the transition out of the Weyl semimetal 
state in the large b limit, i.e. when b = b c 2. Only if 


the condition Eq. (96) is satisfied may we expect to get 


negative magnetoresistance, quadratic in the magnetic 
field in the general case. Otherwise, in magnetic Weyl 
semimetals, linear magnetoresistance, which is allowed 
by symmetry, will dominate the quadratic one at small 
fields. Under the condition Eq. (96), the linear terms in 


magnetoresistance are 0(b/b c 2) and may thus be ignored. 
Linear magnetoresistance is strictly absent by symme¬ 
try, of course, in the case of noncentrosymmetric Weyl 
semimetals. In this case the above results apply without 
restriction. 


V. DISCUSSION AND CONCLUSIONS 

In this paper we have developed a theory of anomaly- 
related weak-field quadratic negative longitudinal mag- 
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netoresistance in Dirac and Weyl metals. An impor¬ 
tant issue is how to differentiate this novel magnetoresis¬ 
tance from other possible contributions, which are more 
mundane in origin. In fact, longitudinal magnetoresis¬ 
tance, which is what we are interested in here, in never 
entirely mundane. The reason is that, from the sim¬ 
plest Drude theory viewpoint, the only possible source 
of magnetoresistance is the Lorentz force, which is of 
course absent for electrons, propagating along the direc¬ 
tion of the field. Drude theory thus predicts that lon¬ 
gitudinal magnetoresistance is always absent, which is 
not the case: there are many examples of materials, ex¬ 
hibiting it, even at low fields. Several possible sources of 
longitudinal magnetoresistance have been identified over 
the years ,®1^] but perhaps the most universal source, re¬ 
lated purely to the intrinsic properties of the electronic 
structure, was described recently by Pal and MaslovP^ 
They have shown that longitudinal weak-field magnetore¬ 
sistance arises necessarily when the shape of the Fermi 
surface exhibits certain types of angular anisotropy with 
respect to the direction of the magnetic field. This 
mechanism gives positive magnetoresistance, increasing 
quadratically with the magnetic field at low fields. For 
the anomaly-related negative magnetoresistance to be 
observable, it needs to be larger than this Fermi surface 
anisotropy-driven magnetoresistance. From this view¬ 
point, the Dirac metal case seems to be the best: one 
may expect both a weak anisotropy and a large axial re¬ 
laxation time in this case. Weyl metal case with either 
a very large separation between the nodes, i.e. separa¬ 
tion approaching the size of the first BZ, or a very small 
separation (unless it arises from a parent zero-gap Dirac 
semimetal), are both problematic, since the axial charge 
relaxation time may be expected to be small in these 
cases. 

It is also important to remember that the theory, pre¬ 
sented in this paper, applies only in the semiclassical 
limit, i.e. ujb/zf ■C 1. In the opposite, ultraquantum 
limit, one may expect a linear negative magnetoresis- 
tancep^^ which may be obtained from Eq. ( |72| by sub¬ 
stituting the lowest LL density of states g{tF) ~ B. It 
is also possible to have quadratic magnetoresistance of 
both signs in this limit, arising from the combined ac¬ 
tion of chiral anomaly and field-induced modification of 
the density of states. 53 One hopes that the current ex- 


perimentf 13 l 22 l 54 t Ml are j n the semiclassical limit, as the 
quadratic magnetoresistance is observed at low fields, but 
the magnitude of the ratio uib/cf is at the moment un¬ 
certain in these experiments. 

Finally, in the theory, developed in the paper, a partic¬ 
ular model of the impurity scattering was assumed: weak 
point-like Gaussian-distributed scatterers. This assump¬ 
tion was made primarily for computational convenience: 
the ladder sum for the diffusion propagator, Eq. (421, 
may only be calculated straightforwardly in this case, 
which is a limitation of the present approach. In many 
cases, a more physically-realistic model should involve 
Coulomb impurities. In this case one might in fact expect 
that the axial relaxation rate should be even smaller, at 
least in the Weyl semimetal case, as the finite-momentum 
scattering, necessary to scatter electrons between the 
nodes, will be suppressed. However, a detailed calcula¬ 
tion of this would certainly be helpful. What happens to 
the magnetoresistance in the strong disorder limiiP^HH 
is also an important and experimentally-relevant issue, 
worthy of a thorough study. 

In conclusion, we have presented a theory of chiral 
anomaly-driven negative quadratic longitudinal magne¬ 
toresistance in Dirac and Weyl metals in the weak mag¬ 
netic field regime. We have demonstrated that this effect 
has two crucial ingredients. One is the coupling between 
the vector and the axial charge density, proportional to 
the magnetic held, which arises from the chiral lowest 
Landau level, or the nontrivial Berry curvature of the 
band eigenstates. This coupling will in principle exist in 
any system with a nonzero Berry curvature and is in this 
sense not specific to Dirac or Weyl metals, although the 
transport coefficient, describing such a coupling, has a 
universal value of eB/2n 2 in Dirac or Weyl metals only. 
The second ingredient is the near-conservation of the ax¬ 
ial charge density, manifesting in large axial charge relax¬ 
ation time T a /To 1. This property is specific to Dirac 
and Weyl metals only and is necessary for the negative 
quadratic magnetoresistance to be observable. 
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